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Abstract 

For the purpose of analyzing non-perturbative dynamics of string theory, Nishimura 
and Sugino have applied an improved mean field approximation (IMFA) to IIB ma- 
trix model. We have extracted the essence of the IMFA and obtained a general 
scheme, the improved Taylor expansion, that can be applied to a wide class of 
series which is not necessarily convergent. This approximation scheme with the 
help of the 2PI free energy enables us to perform higher order calculations. We 
have shown that the value of the free energy is stable at higher orders, which sup- 
ports the validity of the approximation. Moreover, the ratio between the extent of 
"our" space-time and that of the internal space is found to increase rapidly as we 
take the higher orders into account. Our results suggest that the four dimensional 
space-time emerges spontaneously in IIB matrix model. 
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1 Introduction and summary 



After the enormous success in revealing perturbative dynamics of string theory with/without 
D-branes [jl], ^ ^, our focus will go on to the full non-perturbative dynamics, or to the 
theory which has predictive power for the real world. 

In such a stage, both a model which is to be the true theory and a method to extract non- 
perturbative information from it are important. Some models have already been conjectured 
as the constructive definition of string theory and/or quantum gravity ^ i|, 0. Among 

others we are interested in IIB matrix model [^, |T0[ which would be the most promising to 
attack the non-perturbative dynamics of string theory. 



As to the method, J. Nishimura and F. Sugino []rT| have applied an improved mean field 
approximation (IMFA) [0 to IIB matrix model with an excellent idea for breaking the Lorentz 
symmetry, and suggested that the four dimensional space-time is more stable than those of 
the other dimensions. In sec. ^, we review the IMFA with an application to a toy model. 
Although the IMFA is somewhat mysterious, we find that the essential feature of the IMFA 
can be captured in terms of the improved Taylor expansion (ITE) which is elucidated in sec.|[ 
In [ITT] the authors have calculated to the third order and got a result which is both remarkable 



and uneasy: 

remarkable fact Even at the third order, non-perturbative aspects of IIB matrix model can 
be captured definitely. 

uneasy fact Not knowing the stability of the result, we can not strongly believe it. 

To confirm the validity of the IMFA for IIB matrix model, we should investigate the stability 
of the result, especially the free energy, as we increase the order of approximation. The aim 
of this paper is to deepen the remarkable aspect and to remove the uneasy problem at a time. 
So we proceed to the next order calculation. Since Feynman graphs to be calculated increase 
in number and become complicated, we need to introduce a new method. It turns out nice to 
use the two-particle irreducible (2PI) free energy which is explained in sec. |[ With the help 
of the 2PI free energy we have completed the fifth order calculation and see how the above 
aim is achieved. We find that the contribution of the fifth order to the free energy is smaller 
than that of the third order, which supports the validity of this approximation. 

The ratio between the extent of the four dimensional space-time and that of the extra 
internal space is found to increase from three (the third order) to six (the fifth order) as is 
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shown in sec. |^. This indicates that the ratio is infinite in the full theory and the fiat four 
dimensional space-time is generated. 



2 Improved mean field approximation 
2.1 general prescription 

Suppose we have some action function S{x) and its partition function 

Z= /"dxe-^(^). (2.1) 



In general, this integral can not be performed analytically, and we need an approximation 
scheme to evaluate the partition function. In the mean field approximation, we add the mean 
field action S'o(x;a) which can be integrated analytically, and subtract it as an interaction 
term. The mean field action contains a set of parameters a which should be determined later 
according to the order of the approximation. The original action S{x) is also treated as an 
interaction term. Here we introduce a coupling g formally, which will be set to 1 at the end of 
the calculation. To get the k-th order approximation, we convert the expression to a formal 
power series with respect to g and truncate it up to order k. The prescription is as follows: 

4 convert to the formal power series w.r.t. g 
Z r^Yl ^-T~ I ("^(^^ " "^o^^' (2.3) 

n=0 ^' 

J| truncate the summation 
Zk{a) = ^ / da; {S{x) - So{x; a))"e-^°(^^'^). (2.4) 

n=0 ^' 

The truncated partition function Zk{a) is a function of the parameters a, though the 
original partition function is independent of them. Thus we can determine the values of the 
parameters so that the truncated partition function is stationary with respect to a. We call 
this procedure the improved mean field approximation (IMFA). 

To approximate the correlation function, we use the same truncation: 



'^^ > n=0 



■ (2.5) 

neglect 0(g'''+-'^) 



Note that at the first order IMFA is the same as the variational principle. Note also 
that the power series ( p.3|) is in general only an asymptotic series if we fix the parameter 
a. Nevertheless by determining a order by order, we can approximate the original partition 
function with high accuracy as we will see in the next subsection. 



2.2 application to a toy model 

As an example, let us apply the IMFA to the zero dimensional (f)^ model whose partition 
function can be integrated exactly. The action is 

5 = + ^c^'. (2.6) 

The exact values of the partition function and the two point function in the massless case 
{w? = 0, A = 1) in which we are especially interested are 

^masslcss = J d<j) e'' = 2 ■ 3\ ■ (4!)"=^/^ ■ r(l/4), (2.7) 
i^massicss = " log Z ~ -1.389388802, (2.8) 
( )massicss = ^ /^^ ^'e"^ = (Al)'/''^^ ^ 1.655801764. (2.9) 

We adopt the following action as 5*0, 



ml ,2 



^0 = -fr, (2.10) 

where mg is a variational parameter which corresponds to a in the previous subsection. Even 
in the massless case, this action enables us to calculate the free energy = — log Zj^ by using 
the ordinary Feynman diagrams. 

Graphically, Fk is the sum of all connected vacuum graphs shown as in eq. ( 2.11| ), where a 
line and a cross with a circle represent a propagator l/rriQ and a mass counter term g{m'^ — m?) 
respectively. The four-point vertex denotes g\. 




Figs.|I] and |] show the approximated free energies and two point functions as functions 
of ml in the massless case. We find that the higher the order of the approximation, the 
wider the plateau at the exact values in both cases. This phenomenon is expected because 
the free energy F should not depend on the parameter rriQ which is introduced artificially. 
Although the truncated free energy Fk depends on rriQ, we expect ttiq dependence becomes 
weaker for larger values of k. Since the truncated free energies and two point functions are 
merely polynomials in l/rriQ, the existence of the plateau indicates that extremum points of 
these polynomials concentrate on limited regions. Once the plateau exists, small changes of 
ttiq do not affect the resultant values. Therefore if we have a plateau, we can regard the value 
in the region as a good approximation of the exact value. 
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Figure 1: Truncated free energies in the IMFA scheme. The horizon- 
tal axis denotes m^. 



3 Improved Taylor expansion 

The IMFA is a rather curious approximation scheme. It is interesting to understand the 
ground on which it works. In this section, we see that the improved Taylor expansion (ITE) 
explains some of the essential aspects of the IMFA. 
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Figure 2: Truncated two point functions in the IMFA scheme. The 
horizontal axis is mg. 



Suppose we have a sequence of functions {/i(x)}j=o,i,... which does not necessarily converge. 
We define an improved sequence xo)}i=o,i,... by 



fr-^'ix- Xo) = j2-^^- ^oTtiM- (3.1) 



n=0 



We call this prescription the improved Taylor expansion (ITE). We emphasize that each 
ymprovcd constructed from finite number of //s, that is, /o, /i, . . . , /«. This differs from 
the ordinary Taylor expansion in that we use fj with smaller j's for higher derivative terms. 
As we see below, different values of Xq make {/i™^'^°™*^}j=o,i,... converge in different regions. 
This property enables us to obtain a series which converges around various values of x by 
choosing proper xq correspondingly. 

For example, we consider a power series given by 



- , X 

n=0 



and define {/i}i=o,i,... by 



n=0 
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The behavior of this sequence {/i}i=o,i,... is shown in Fig.|^, which manifests that the conver- 
gence radius of the series is 1. The improved sequence a;o)}j=o,i,... for xo = 1.1, 
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Figure 3: Original series. The horizontal axis is x. 



for instance, behaves as in Fig.^. This sequence converges well around x ~ 3/2. Thus the im- 
proved sequence has a totally different convergence region from that of the original sequence. 
As we can see from this example, the reason why we use fj with smaller j's in the definition of 
the ITE eq. (|3.1| ) is that fj with the higher j's make behavior of the sum worse for a divergent 
series. 

To see the relation with the IMFA more clearly, let us fix x to 3/2, outside the convergent 
radius {x corresponds to in eq. (p.6|)). Instead we regard /™''™™*^(3/2; Xq) cl-S Sb function of 
Xq. Now Xo corresponds to a parameter in the IMFA {rriQ in eq. ( p.lO|) ). As seen from Fig.^, 
there are plateaus on the exact value /(3/2) = 0.4. The plateau gets wider as we increase the 
order of the approximation. The approximation is nice and the improved functions are almost 
independent of the variational parameter xq as expected. 

The expansion ( p.l|) can be written in a concise form. First, we introduce a formal expan- 
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Figure 4: ITE around Xq = 1.1. The horizontal axis is x. 
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Figure 5: The plateau in ITE scheme. The horizontal axis is Xq. 
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sion parameter g, and define a function from the sequence of functions {/j(x)}i=o,i,. -' 

oo 

f{x) = Y,9'H^), (3.4) 



i=0 



where hi{x) = fi{x) — (z = 1,2, . . .) and ho{x) = fo{x). For example, in the case of 

power series eq. ( |3.2| ) we have hi{x) = (— )* x\ Next we substitute xq + g{x — xq) for x in this 
f{x) and truncate it with respect to g: 

i 

Here |j denotes neglecting 0{g''~^^) terms and setting g=l. One can easily verify that this 
expression reproduces the functions in the improved sequence eq. (|3.1|). 

In relation to the IMFA, the derivatives correspond to insertions of mass counter terms. 
As a result, the IMFA is included in a more general scheme, the ITE. The abstract feature of 
the ITE may be useful to find the essence of the IMFA. Once we understand the essence, we 
could refine the IMFA in various ways and use a refined method to grasp the non-perturbative 
properties of IIB matrix model further. 

For reference, some curious behavior of the ITE and an application to an asymptotic series 
are examined in appendices and 0. 



4 2PI free energy and Legendre transformation 

In the application of the IMFA to various models, we can drastically simplify the computation 
by using the 2PI free energy. In this paper, we apply the Legendre transformation to the 2PI 
free energy to obtain the ordinary free energy, and search its extrema. 

Although within this procedure the 2PI free energy seems just a technical tool, we can give 
it a more fundamental meaning related to the Schwinger-Dyson equations (SDE). In fact, we 
can apply the ITE scheme directly to the 2PI free energy, which we explain in appendix 0. 



4.1 Schwinger-Dyson equations and 2PI free energy 

The SDE are the consistency conditions among the Green functions and in some cases we can 
extract non-perturbative information from them. Again we consider the massless 0^ model. 
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For the full propagator c = ( 0^ ), we can derive the SDE as shown in eq. ([4.1|) . 




The derivation of this equation is given in appendix 0. Note that each graph has no self-energy 
part when we regard each full propagator as a single line. 

This equation is, in fact, the stationary condition of the two particle irreducible (2PI) free 
energy which we call G hereafter. 




2PI graphs are the ones that remain connected after cutting any two distinct propagators. In 
other words, 2PI graphs do not contain any self-energy part as a subgraph. This property 
comes from the fact that the SDE do not contain any self-energy part. The correspondence 
between the SDE for the full propagators and the 2PI free energy holds generally in field 
theory. 

The advantage of considering the 2PI free energy is that we can greatly reduce the number 
of graphs. Next, we will explain the relation between the 2PI free energy and the ordinary 
free energy. 



4.2 Legendre transformation 

Let us start with the free energy, 

Z{J) = [d(f)e-^^'^^-^'^\ (4.3) 



F(J) = -logZ(J). (4.4) 

Here we introduce a source J for 0^. We apply the Legendre transformation to F{J) with 
respect to J. 

KiJ) ^ (4.5, 
9 



G{K) 



F{J) - J 



dFjJ) 



(4.6) 



J=J{K) 

where K is the conjugate variable of J. To investigate the original theory, we set J = 0. Then 
K{J = 0) is equal to the full two-point function c, 



K{J = 0) = {4>')j=o = c, 



(4.7) 



and G{K = c) is nothing but the sum of the 2PI vacuum graphs |]T3 



For example, we use this procedure for the 0^ model, eq. (p.6|) . The free energy F is given 
by the sum of the connected vacuum graphs, 

A 1 



log(m^ + 2 J) 



2jy 



1 1 

48 ^ 16 



1 



2jy 



0{X' 



(4.8) 



where each term corresponds to the one in eq. (|2.11|) other than the graphs with mass counter 
terms. The full propagator becomes as follows: 



1 



A 



+ A^ 



+ 0(A^ 



(4.9) 



2(m2)3 3(m2)5 

This equation can be solved with respect to 1 / iteratively, and we get the 2PI free energy 
G. 



1 



G 



c + -r H 

2 12 

1 



1 o A 



logc H — cm" 

2 2^ 2 



A^ 
48 



(4.10) 
(4.11) 



which reproduces eq. ( |4.2|) correctly when = Q. 

In the practical calculations, it is easier to compute the 2PI free energy G{c) first and then 
obtain F{m?) by applying the inverse Legendre transformation to G{c). 

Once we obtain the free energy F{m?), we replace A with g\ and apply the ITE with 
respect to im? . For example if we consider the massless case corresponding to IIB matrix 
model, the improved free energies are 



limprovcd 



F{ml-gml)\^ 



log CQ + g 



A 



A2 
48^ 



A^ 



16"° 



A 



+ 



(4.12) 
(4.13) 



where cq stands for l/mg. Again, the last equation is the same as eq. ( p.ll|) term by term. 



"'^Here an additive constant is not important. 
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5 The ansatz and the results 



5.1 IIB matrix model 

Our aim is to analyze IIB matrix model [^]. The action is 

S = -^Tr( \[A„A^r -li^^n^^^P] ), (5.1) 

where = 1, . . . , 10) and ip°'{a = 1, . . . , 16) are all x hermitian matrices transforming 

as a vector and a left-handed spinor representation under SO(IO). The symmetries of IIB 
matrix model are the matrix rotation U(A^), ten dimensional Lorentz symmetry SO(IO) and 
the type IIB supersymmetry. As discussed in |TD[ , the distribution of the eigenvalues of is 
interpreted as the space-time itself. 

It is difficult to integrate this action exactly. So we apply the IMFA by adding and 
subtracting a quadratic term. We specify the full propagators for bosons and fermions instead 
of the quadratic term Sq. Now we focus on the dimensionality of space-time or the ratio 
between "our" space and the "internal" one, and consider such ansatz that the U(A^) symmetry 
remains unbroken. So gauge indices are set to preserve the U(A^) symmetry. 

{A^'^A^\)=g',C^,^)6\S^', (5.2) 
( / ) = g'o ^n,.p]{r'''C-'r^ ^ (5.3) 
where C is the charge conjugation matrix such that 

*C = -C, CTf" = -'Tf'C {fi = l,...,W). (5.4) 

Since the U(A^) index part is symmetric in each propagator, the Lorentz index part is forced 
to be symmetric for bosons and anti-symmetric for fermions. We have converted the spinor 
indices a, (3 to the rank three anti-symmetric tensor uy^^p] through the gamma matrices, which 
is equivalent to the antisymmetric product of two left-handed spinor representations. 

Here we give an overview. First, we calculate the 2PI free energy G by evaluating planar 
vacuum diagrams with the above exact propagators. We consider the large- A^ limit with the 
't Hooft coupling g = g^N fixed to 1. We use this g as the formal expansion parameter in the 
ITE later. Secondly, in order to investigate in what shape the eigenvalues distribute and to 
compare the free energies of various vacua we try various ansatz with symmetries which are 
subgroups of SO (10). With these symmetries we can restrict the form of the propagators and 
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reduce the number of variational parameters. Thirdly, for each ansatz we make the inverse 
Legendre transformation to the 2PI free energy G to get the ordinary free energy F. We 
further apply the ITE with respect to the conjugate parameters of C and u as in eq. ( [4.12| ) 
and obtain the improved free energy 

^improved _ Finally, in each ansatz we extremize F^^p™™'^ 
and compute the extent of the space-time from ( Tr (A^j^A^) ). 

5.2 2PI free energy 

We calculate the 2PI free energy G to the fifth order. The Feynman rules and all graphs 
with their symmetry factors are shown in appendices and 0. At a glance we notice that 
the number of graphs are much reduced by the condition of 2PI. For example, at the third 
order there are only four 2PI graphs as in appendix |E.3| , although thirty four connected graphs 
contribute to the ordinary free energy. 

The general form of the 2PI free energy is given hyf\ 

= 3(l + log2)-ltrlogC + ^trlog^-| (trC^ - (trC)2+tr(^r'^^r-)C^,) + ■ ■ ■ , (5.5) 

where -^"^ = u^,yp(r'^'^'')"^/3! and the traces are taken in the vector representation for G's and 
in the left-handed spinor representation for i^'s. Although we can write down the higher order 
terms without evaluating the traces as in eq. ( p.5| ), it is technically difficult to evaluate these 
traces and obtain a concise expression of G. Instead we evaluate the traces for each ansatz. 
We present in appendix |G| the explicit forms of G to the fifth order for the SO (7) and the 
SO (4) ansatz which are explained in the next subsection. 

5.3 symmetries and ansatz 

The guideline of making the ansatz is as follows: 

• First of all, we fix a subgroup of SO(IO) from S0(7) to S0(1), which will be interpreted 
as the rotational symmetry of "our" space-time. The component of the fermion two- 
point function u is non-zero only when no indices are along "our" space-time by this 
symmetry. In fact, we shall see that the extent of "our" space is larger than that of the 
remaining "internal" space. 
^We adjust an additive constant to the definition of 



12 



• Because m is a rank three anti-symmetric tensor, we decompose the internal space into the 
sum of three dimensional spaces and the rest, e.g. 4=3+1 or 8=3+3+2. We then assume 
the rotational symmetry for each space, e.g. S0(3)xS0(l) or SO(3)xSO(3)xSO(2). 
This symmetry forces the component of u to be zero unless all three indices lie in one 
of the three dimensional spaces. 

• Finally, we impose the permutation symmetry on the internal three dimensional spaces. 

In this procedure, it is very important to include all parameters permitted by the imposed 
symmetry. Even if the free energy has an extremum in a randomly restricted parameter space, 
it need not be a real one in the full parameter space. On the other hand if we restrict the 
parameters by a symmetry, the free energy is promised to be stationary in the directions of 
the parameters which break the symmetry. 

SO (7) ansatz 

The fermion propagator is represented by the rank three anti-symmetric tensor u. There- 
fore if we give a non-zero value to only one element of u, SO (10) symmetry breaks down to 
SO(7)xSO(3). Under this symmetry group, propagators are represented by three parameters, 

Vl,V2,U. 
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/ V, 






\ 


7 
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10 


\ 






V2 ) 



i = ^r^'S'io. (5.7) 

so (6) ansatz 

Next, we consider the case where S0(6) symmetry survives, uy^^p^ has non-zero values only 
in the extra four dimensional space. A tensor u^^pp of SO (4) is dual to a vector u^. Therefore 
up to SO (4) rotation, all configurations are identical. As a result SO (3) C SO (4) transverse to 
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remains unbroken. Therefore this ansatz has SO(6)xSO(3) symmetry. 
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( 
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Cfj^u = 7 
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V, ) 



(5.8) 



^ = ur«'9'io. (5.9) 

SO (5) ansatz 

With five extra dimensions, the situation becomes a httle comphcated. u^^yp\ of SO (5) is 
now dual to a rank two anti-symmetric tensor (adjoint representation) u^^'^y In general, u^'^'^^ 
can be taken to the standard form by SO (5), there exists two free parameters. For simplicity 
(though it may not be natural), we set one of them to zero. As to the fermion propagator, 
after all, we take the same form in these three cases. This ansatz has S0(5) xS0(2) xS0(3) 
symmetry. 

C^u = diag(five Vi's, two V2's, three Vs's), (5.10) 
^ = (5.11) 



SO (4) ansatz 

Assuming S0(4) symmetry, the extra dimension is 6 = 3+3. We can naturally impose 
simple symmetry S0(3) xS0(3) xZ'2 where Z2 symmetry exchanges two S0(3) factors with 
the reversion of the 1st direction (for no parity change in total). Consequently we arrive at 
the three parameter ansatz. This ansatz has S0(4) xS0(3) xS0(3) XZ2 symmetry. 
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( 
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V2 
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V2 ) 



V2 



(5.12) 



(5.13) 



SO (3) ansatz 
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We decompose the extra 7 dimensions to 3+3+1 , and impose S0(3) x S0(3) x symmetry. 
The fermion propagator is the same as the above case. This ansatz has S0(3) xS0(3) xS0(3) XZ2 
symmetry. 

C^jy = diag(three Vi's, one V2, six Vs's), (5-14) 

^ = _L„(r5.6,7 + r«,9,io)_ (5 15) 

In the first S0(3) factor (which will be interpreted as "our" space-time), the reversion of 1st 
direction included in Z2 forbids the non-zero value of 1*1,2,3- 

SO (2) ansatz 

We decompose the extra 8 dimensions to 2+3+3 , and impose S0(2) xS0(2) xS0(3) xS0(3) x 
symmetry. We then have the following form: 

C^y = diag(two Vi's, two V2S, six Vs's), (5.16) 
v2 

S0(1) ansat^ 

Finally we consider the one dimensional space-time. We divide nine extra directions into 
three parts each of which has three dimensions. Three SO (3) symmetries and permutation 
of these factors are imposed. S3 is combined with the reversion of the 1st direction in order to 
be a subgroup of the original SO (10). This ansatz has S0(3) xS0(3) xS0(3) x^a symmetry. 

C^iy = diag(one Vi, nine V2's), (5.18) 
V3 

5.4 free energy 

The ansatz listed above have three or four parameters. Namely the boson propagator has two 
or three parameters Vi and the fermion propagator has one parameter u. After substituting 
each ansatz into the general expression of G as eq. ( [5.5| ), we obtain the 2PI free energy G{Vi, u). 
The explicit forms of GiVi, u) are given in appendix [G] for the S0(7) and the S0(4) ansatz. As 
we discuss later, these cases turn out to be of particular importance. Then we make the inverse 
■^We use the term S0(1) to indicate one dimensional space-time. 
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Legendre transformation to G{Vi, u) in the following way. We define conjugate variables to Vi 
and u as 

M^ = -^G{y,u), (5.20) 

m = -^G{y,u). (5.21) 
ou 

The ordinary free energy F is 

F{M\ m) = G - ^ M% - mu. (5.22) 

i 

Then we apply the ITE to F{M\m): 

^mproved^^.^ m; M^, mo) = F(M^ + g{M' - Mi,), mo + g{m - mo)) . (5.23) 

k 

Here \k denotes expanding with respect to g, neglecting 0{g^'^^) terms and setting g = 1 a.s 
in eq. ( |3.5| ). Since the action of IIB matrix model does not have a mass term, the improved 
free energy is given by 

^improved^Q^ 0; M^, mo) = F(M^ - gMi,, mo - gmo) . (5.24) 

k 

The explicit forms of F{Mq — gMQ^rriQ — grrio) for the S0(7) and the S0(4) ansatz are given 
in appendix |^ . We differentiate this free energy with respect to each parameter and set the 
results to zero. By solving these equations we get the extrema of the improved free energy. 
If there are more than one solutions, we select the one which gives the smallest free energy 
among them. Such a solution is expected to be on the plateau. 

Fig.|^ shows the free energies for the SO (7) and the SO (4) ansatz. As we explain below, the 
other cases are reduced to these two cases. Their numerical values are shown in Table |l|. The 
corresponding values of Mq and mo are shown in appendix |. We find that the contribution of 
the fifth order is smaller than that of the third order. This fact suggests that the IMFA works 
well, and we can trust these results as in the case of the 0^ toy model. 

Some comments are in order. 

• At the 2nd and the 4th order there is no reasonable extremum. As in the case of the 0^ 
model, extrema do not necessarily appear in the even and lower order levels. Therefore 
we consider only the odd orders. 
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ansatz 


quantities 


1st order 


3rd order 


5th order 


S0(7) 


F 


5.522722156 


1.620942535 


1.492429842 


R 


0.6407423748 


0.7670244718 


0.8051803940 


r 


0.3276944430 


0.3556261073 


0.3844627899 


P 


1.955304364 


2.156828390 


2.094299930 


S0(4) 


F 


6.153347114 


0.696885300 


0.077693647 


R 


0.7500535231 


1.161504610 


1.7090706440 


r 


0.4038451083 


0.3796625830 


0.2691460598 


P 


1.857280202 


3.059307619 


6.349974602 



Table 1: Numerical values. 
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• For the S0(6) ansatz, the extremum solution reduces to that of the S0(7) ansatz, i.e. Vi 
and V2 in eq. (|5.8| ) take the same value and the SO (6) symmetry is enhanced to SO (7). 
In the same way the SO (5) ansatz reduces to the SO (7) ansatz. This is natural because 
the fermion propagators have the same form in these cases. On the other hand the SO (3) 
and the SO (2) ansatz reduce to the SO (4) ansatz. We thus conclude that the fermions 
play a crucial role in the process of the spontaneous breakdown of SO(IO). 

• The S0(1) ansatz has no solution even at the 1st and the 3rd order. This may indicate 
that one dimensional space is not realized as "our" space. 



5.5 extent of space-time 

Next, let us consider the extent of "our" space-time R and that of the "internal" space r. 
Here, for both the SO (7) and the SO (4) ansatz, R and r are defined by 

9 ,1 , o , dF(M\m) 

R'={^ TrA,^ } = . (5.25) 

9 ,1 . 9 ^ dF(M\m) 

= { ^ TVA,„^ > = L_^. (5.26) 



These quantities can be also improved by the ITE, and obtained as in eq. ( ^.24|) : 



1 (9F 

( - Tr V )k = -Qj^iK - gM'o, mo - gmo) ^, (5.27) 

1 dF 
{ - Tr V )k = -QM^iK - 9K, ^0 - gmo) ^. (5.28) 

We evaluate these quantities using the values of Mg and rriQ given in appendix |. R and r are 
shown in Fig.|^ and their ratio p = R/r is plotted in Fig.|^. 

In contrast to the stability of the free energy, both R and r for the SO (4) ansatz changes 
considerably. It is striking that their ratio grows rapidly as we increase the order of the 
approximation. 

The fact that the SO (4) ansatz gives the lowest free energy and the large value of the 
extent ratio indicates that the path integral of IIB matrix model is dominated by the nearly 
four dimensional distribution of the eigenvalues. 



6 Conclusions and discussions 



Our main results are the following: 
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S0(7) r 






S0(4) R * 
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S0(4) r B 
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1st order 3rd order 5th order 

Figure 7: Extent of the space-time. 
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S0(7) 
S0(4) 




1st order 



3rd order 



5th order 



Figure 8: Ratio between the extent of "our" space-time and that of the "in- 
ternal" space. 
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• Four dimensional configurations dominate in the path integral of JIB matrix model. 

• The fermionic field plays an important role in the breakdown of the Lorentz symmetry. 
Its two point function seems to determine the remaining symmetry. 

• For the SO (7) ansatz, the free energy and the extent ratio become stable at the 5th 
order. On the other hand, for the SO (4) ansatz, the free energy still decreases at the 5th 
order, but it seems to approach to a definite value. On the other hand the ratio grows 
more rapidly as expected. Stability of the result requires still higher order analysis. 

The IMFA is, in a rigorous sense, not an "approximation" scheme. We do not know the 
true reason why a free energy converges to its exact value. But if we apply the IMFA to 
some models, the results are magically nice. To find a trick behind the IMFA should be very 
important. 

In any ansatz except the S0(1), the "internal" space shrinks in such directions that the 
fermion full propagator has non-zero value. We need to understand this phenomenon, shrinking 
by fermion, at least qualitatively. From this point of view, it would be interesting to clarify 
the relation between our result and the mechanism of the symmetry breaking by the phase of 



the fermion determinant suggested in |T5 



In relation to the above connection between the "compactification" of the space and the 
fermion full propagator, we make a brief comment on a possibility that the full SO(IO) sym- 
metry will be restored. One might wonder that eigenvalue distribution having the full SO(IO) 
symmetry is likely to be realized. However, as far as our approximation scheme is concerned, 
such distribution corresponds to setting all of the components of u to be zero, in which both 
the ordinary free energy and the 2PI free energy apparently diverge logarithmically. Therefore, 
it cannot be a minimum energy configuration or a solution to the SDK. This is consistent with 



the results in 15 



Growth of the extent ratio for the SO (4) ansatz is quite interesting. To what value will it 
reach as going to the higher order, infinite or finite? If it is approaching infinity, we can regard 
it as an indication that the space-time is spontaneously compactified to the four dimensional 
flat space. 

At any rate, the IMFA partly grasps the non-perturbative dynamics of string theory. We 
have extracted information on the dimension of the space-time as the first step. Probably 
the next step is to predict the gauge symmetry. Although some refinement to the IMFA is 
required, we will become able to achieve it in the near future. 
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There may be some directions to refine the IMFA. 

• Not only two point functions but three or more point functions can be included to the 
SDE. Also it may be interesting to include one-point functions. 

• In the present paper, we have used the 2PI free energy and the full propagators just 
as a calculational technique to reduce the number of graphs. But the notion of the full 
propagators, which has a direct connection with the SDE, may be more powerful in 
examining various aspects of the non-perturbative dynamics. In light of this, it would 
be important to investigate a direct application of the ITE to the 2PI free energy as 
done in appendix |C[ 

Although there remain many unknown or ambiguous problems such as the existence of 
gravity, chiral fermions, the connection with the low energy effective field theory and so on, 
we hope we can proceed step by step with this new method. 
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A Some curious behavior of ITE 

In sec.^ we apply the ITE to the series "^^^ value at a; = 3/2 is well approximated. 

The method in sec.^ however, is not effective when we are to evaluate the series farther from 
X = 1. Therefore we use the ITE in a different way. The function we consider is 




(A.l) 



fk{x) = f{x)\^. 



(A.2) 
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As X becomes large, the behavior of the power series in x would be worse. In such a case, a 
more proper variable is the inverse of x, and we apply the ITE to it, i.e. 



1 1 

— + 

X Xq 



f, 



improved 



/ 



1 1 

X Xq 

1 



(A.3) 
(A.4) 



The original series of j[x) without the ITE behave as shown in Fig.|^. 

After applying the ITE, the behavior changes drastically. Fig.|^ shows the improved series 
at Xo = 0.8. The improved series behaves well in the region a; > 1 in which the origi- 





1 1 1 

order 5 




order 10 




\ order 15 




i order 20 




!\ order 25 




|\ order 30 


/ 
/ 


i / ~~ ^ 

1/ — ~ 

/ 
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Figure 9: Improved series. The horizontal axis is x. 



nal series diverges. For example, if we evaluate the value at x = 3, the error| is almost 

■j^g — (order of approximation) 



. . , n 1 , (approximated value) — (exact value) 

horror is denned here as 



(exact value) 
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B Improved Taylor expansion for asymptotic series 



Although our actual interest is in analyzing IIB matrix model and its planar sum has a finite 
convergence radius, ITE can also be applied to asymptotic series. We have already seen in 
sec.0 that the IMFA for the (f)^ model works well. The expansion in g of this model is only 
asymptotic and not convergent. As another example, we apply the ITE to F function here. 



J{x) = ( logr(x-^/2) -x-^/2(logx-i/' - 1) + - log 



X 



-1/2 



2lT 



1 

12 



360^^^ 1260 



2 2 

9 X 



4 4 

;9 X 



691 



-g^x^ 



(B.l) 
(B.2) 



1680" 1188" 360360^ 
which is also an asymptotic series. As in the above example, we apply the ITE for the inverse 



variable. Figs.|I3 and |11| show respectively the truncated functions without improvement and 
with the ITE around xq = 0.03. The behavior in the large x region is smoothed out after 



0.085 



0.084 



0.083 



0.082 



0.081 



0.08 



exact function 
order 5 
order 10 
order 15 
order 20 
order 25 
order 30 




0.2 



0.4 



0.6 



0.8 



Figure 10: Original truncated functions. The horizontal axis is x. 



the ITE. If we fix x to 1/3, conversely, there is a plateau for each improved function as shown 
in Fig.[T2|. The errors from the exact values at the extremum are about 10^^ for the improved 



functions of order between 6 and 26. On the other hand the truncated functions without 
improvement have errors of about 10~^, if the order of truncation lies between 3 and 7. As 



23 



0.085 



0.084 



0.083 



0.082 



0.081 



0.08 




0.2 



0.4 



0.6 



0.8 



Figure 11: Improved functions. The horizontal axis is x. 
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Figure 12: Plateau for each improved function. The horizontal axis is 
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going to larger values of x, errors grow in both cases. At any rate the error becomes small 
and the oscillation becomes smooth by the ITE. 



C Direct application of ITE to 2PI free energy 

The ITE can be applied not only to the ordinary free energy but also to the 2PI free energy 
which has a direct connection to the SDE. 



Let us consider the massless matrix 0^ model [jT4 



^ = iV^^Tr0^ (C.l) 

In the planar limit the perturbative expansion has a finite convergence radius. Therefore we 
expect that this model is more tractable than the usual (p^ model. 

The 2PI free energy divided by A^^ is expanded as follows : 

G{x) = -\ log(x) + - + ^a:^ - + 1^ + 0(/), (C.2) 

where x is the full propagator. We can improve this function by replacing x with XQ+g{x — XQ). 
If we fix X and regard G as a function of Xq, we see a clear plateau as shown in Figjl 
When we plot the improved 2PI free energy G by two variables x and xq as in Fig.p 



the contour at G = 0.4 runs in parallel with xq axis at x ~ 0.77. Consequently, we have the 
plateau for xq and achieve a nice approximation around there. Here, in contrast to the case of 
the ordinary free energy, it is the full propagator x which gives the minimum of G that we are 
interested in, because it is nothing but a solution to the SDE. Fixing xq = 0.585 and plotting 
G as a function of x, we succeed to evaluate x = 0.77 at the minimum of G, Fig.p!5|. This 
result almost reproduces the exact value 4v^/9 ~ 0.7698. In this way we are able to apply 
the ITE directly to the 2PI free energy. There is no need for Legendre transformation, and 
we have more direct connection with the SDE. 

D Feynman rules 

We write down the Feynman rule for IIB matrix model. 
• Propagators: 
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0.2 0.3 0.4 0.5 0.6 0.7 0.8 



Figure 13: Plateau of improved G for x = 0.8. The horizontal axis is Xq. 



order 1 2 
0.6 




0.8 



Figure 14: 3D plot of improved G at the 12th order. (0.2 <x< 1.2, 0.8 < 
xo < 1.2). 
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order 15 
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Figure 15: Minimum of improved G. The horizontal axis is x. 
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Extra factors: 

— —1 for each fermion loop. 

— N for each color loop. 



E Planar Feynman graphs 

Here we show all the Feynman graphs from the zeroth order to the fifth order. The number 
below each graph is the symmetry factor. The value of each graph is given by 

(—1) X X (factor from Feynman rules). (E.l) 

symmetry factor 

The first factor of minus one originates from the minus sign of the definition of the free energy, 
F = - log Z. 

E.O zeroth order 




E.l first order 




E.2 second order 




8 4 
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E.3 third order 




29 



E.5 fifth order 
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1 2 
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1 2 



F Derivation of SDE for full propagator 

To begin with, we write the full propagator by the IPI propagator. 




By multiplying the inverse full propagator from the right and the inverse bare propagator from 
the left, we get the equation, 

(• •)"' = ( + -@— (F-3) 
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After combining eg. ( [F.3| ) with the relation of IPI propagator with the full propagator, 




we are able to take the massless limit (set the inverse bare propagator to zero). As a result 
we get the eq. ( [4.1D . 



G 2PI free energy 

Here we present the 2PI free energies for the SO (7) and the SO (4) ansatz. 

^ = + 3 + 3 ln(2) - 1 ln(y7 V^) + \ Hu^') 

+g (21 + ?>V^ + 21 Vi V2 - 56 Vi + 24 V2) 

- ^ - f Vl - 140 Vl - 168 V^i - 12 V^) 
+^^(126 Vl + 14 + 63 + 21 + 91 - 672 

- 96 - 1400 - 168 Vi - 504 V2 - 232 V^) 
+^(-672 vi - 168 vi Vl - 756 V^ V^ - 210 Vl V^ - 630 Vl - 1134 

- 75 vi - 840 vi vi - 168 Vi Vi - 336 Vi Vi + 5040 Vi 

+ 528 1/2*^ - 8064 V^ + 2688 Vi + 9744 Vi Vi + 4032 Vi Vi 

- 22372 Vi - 22512 Vi V2 - 4032 Vi Vi + 1344 Vi Vi 

+ 168 n« Vi Vi - 12096 V2 Vi + 924 u« Vi) 

+g^{l2mQVi^ViVi 

70686 

+ 8316 Vi Vi + 1^/° + 9072 Vi Vi - 24672 Vi - 63504 Vi 

+ 10206 ^8 Vi + ^ + V^i^ + 17856 u^^ Vi - 392000 Vi 

+ 22176 M*^ vi - 4080 1/2^ + 1764 Vi Vi + 1470 Vi Vi + 3696 Vi Vi 
+ 3072 vi - 96768 - 5712 Vi Vi - 498624 

- 279552 1^2^ - 139776 u^^ Vi Vi - 22176 Vi Vi 

+ 23184 Vi Vi - 9744 Vi Vi - 24528 Vi Vi + 27216 Vi Vi 
+ 90720 V2 Vi - 96768 Vi Vi + 2016 Vi Vi + 11424 Vi Vi 
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+ 672 Vl - 3024 - 64512 V{' V2 - 134400 

- 86016 - 59136 V^), (G.l) 

= -- \n{V^ V^) + 3 + 7 ln(2) + - ln(— u^^) 
m 2 ^ ^ ^' ^ ' 2 ^256 ' 

+ g{<o Vl + 15 Vl + 24 Vx V2 - 32 Vi) 
45 

+ (-9 - y ^2" - 36 - 32 - 192 Vi V2 + 96 V^) 

(30 + 85 + 36 l^/ + 60 + 104 - 192 

- 384 V^i^ ^2 ^ ^4 ^4 _ - 768 m*^ V^^ V2 - 1152 V^i 
-256 u^lf ) 

+^(-660 - 324 \/2^ - 600 - 360 

- 189 - 690 If - 822 + 1152 I^^ + 768 1// Vf 
+ 1536 Vi - 768 Vf - 1152 V{' - 4608 V2 

- 3840 Vi - 3072 Vi Vl + 384 Vi Vi + 2304 u"^ V^ 

- 8448 Vi V2 - 8448 Vl 1^2^-16128 Vi V^ - 1792 V^ 
+ ^M<du^Vi) 

+/(-121344m1° V, Vi + 3456 V^ Vi + ^ F^^^ 

+ 6264 Vi Vi - 26880 Vi - 10368 Vi + 3240 Vf Vi + 7860 ^2^° 
44^44 

+ Vi Vi - 43008 - 15872 + 1152 Vi + 7296 

+ 7200 Vi Vi + 8160 Vi Vi + 9480 Vi Vi - 84480 Vi + 4608 V^^^ 

- 12288 Vi Vi - 122880 u^^ Vi V2 - 245760 V^^^^ 

- 153600 l/i^ ^3 _ 7296 v!^ Vi Vi - 12288 Vi Vi + 11520 Vi Vi 
+ 20736 vi vi + 27648 V2 Vi + 102912 Vi Vi - 9216 Vi Vi 

- 12288 Vi Vi + 3072 Vi Vi - 7680 Vi Vi - 32256 Vi V2 

- 1 13664 Vi Vi - 92160 Vi Vi - 49152 Vi Vi). (G.2) 



H Improved free energy 

We present F(M^ - gM^mQ - gnio) for the S0(7) and the S0(4) ansatz. The F™™"^ 
for each ansatz is obtained by neglecting 0{g^^^) terms and setting g = 1 m the following 
expression. Here we use variables Pq and go for notational simplicity, which are defined by 
= -7/(2Moi), P2 = -3/(2M2) and go = 8/mo for the S0(7) ansatz and = 
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Pq = — 3/Mq and go = 8/mo for the SO (4) ansatz. They represent the propagators of the 
mean field action 5*0. 

^Fso{7)(M^ - gM'o, mo - gnio) = 

-iln((Po^)^(Po^f) + 3 1n(2) + iln(go^6) 

+ 9{2Aql (P^) + 3 (Po^)^ + 21 (P^)' - bQql (Pj) + 21 {P^) [P]) + 3) 
+/(-^ {P'of {P'o? + 42 (Po^) (Po^) - f (Po^)^ 

- 96 ql {P^f + 336 ql {P^f (P^) - 84 (P^f (P^) - 168 ql (P,') + 42 (P^)' + 6 (p2)2 + ^ 

^67 

- 60 {P^r + 196 go' (Po^)^ - ^ (Po^)^ - 252 {P^f (Po') + 72 go^ (Po^) + 672 q^ {P^f 

- 336 go^ (Po^)^ (Po^) - 840 go^ (Po^) (Po^)) 

+ /(I - 18144 g2 {P,')' + 1680 go^ (P,')' (P^) + 6615 (P^)' {P^f + 4032 ql [P^f [P]) 

- 5040 go^ (Po^) (Po^) - 12096 go^ (Pj)^ (Po^) + 11424 {P]f {P',f q^ + 5040 go^ (Po^)^ (Po^) 
+ 15624 g6 [P^,f (p2) - 1680 go^ [P^f (P^) + 6384 g^ (P^f (P^f - 3Q9Qql (P^f {P^f 

- 168 go^ (Po^) (Po^)^ + 1152 qt {P.^ + 7728 go^ (P^f (P^) + 2541 {P^Y {P^Y + 3976 ql {P^f 
+ 3360 go' {P^f + 5747 (Po^)^ (Po^)^ - 336 go' (Po^) + 1176 go' (Po^)' + 6804 (Po^)^ (Po') 

+ 63 (Po') (Po^) - 966 {P^r {P^r + 1056 ql {P^f - 336 {P.^ (P^) + 144 go' (Po') - 480 ql {P^f 

- 360 go^ (Po')' - 1256 go^ {P^f + 924 (P^) {P^f - 1008 {P^)' (P^) + 166 {P^f + 6678 {P^)' 
+ 9 (Po')' - 66 (Po')^ - 1134 (Po^)^ + 63 (Po^)') 

+^7^(64176 ql {P^r {P^f - 127008 ql (P^)' + 314496 ql {P^f {P^f 

- 313488 go^ {P',Y (P^) + 5040 ql (P.^ (P^) + 39690 {P.^ (P^)' + 28224 ql (P^)' (Po') 

- 17640 go^ (Po^) (Po') - 255528 (Pj)^ (Po')' go' - 182448 go^ (Po^)^ (Po')' 

- 191688 go^ (Po')^ (Po^)' - 84672 ql (Pj)^ (P^) + 91392 (Pq^)' (Pq')' go^ + 40320 go' (Po')^ (Po^) 

- 194544 go' (Po^)^ (Po')^ + 140616 go^ (Po^)' (Po') + 108360 g« (Po^)' (Po')' - 96642 (P^)' (Pj)^ 

- 272160 (Po') (Po^)^ go' - 5040 go' (Po')' (Po^) + 544320 go' (Po^)^ (Po') - 100212 g« (Po^)^ 
+ 9744 go^ (Po^) (Po')^ + 44688 go' (Po')^ (Po^)' - 25872 go' (Po^)^ (Po')' - 1512 go^ (Po^) (Po')' 

- 13104 g« {P^f (Po') + 169680 g« (Po^) (Po')^ + 9216 go^ (Po')^ - 73920 go' (Po')^ (Po^) 

+ 61824 go^ (Po^)^ (Po') - 101472 go^ (Po')^ (Po^) - 261072 go^ {P^f {P^ f - 231840 go^ (Po^)^ 
+ 14220 ql {P^f - 109200 go^ (Po')^ (Po^)' - 147840 go' (Po')' (Po^)' + 15246 (Po')^ (Po^)' 
+ 35784 g^ (Po^)=^ + 10080 ql {P^f + 34482 (Po^)^ (Po')' - 560 ql {P^) - 26448 g^ (Po')' 
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+ 4116 {P',f - 187110 {P',Y {P^f + 40824 {P^,f (P^) + 84 {P^) (P^) + 13920 (p2)5 

- 2415 {P^f {P^f + 7392 q', {P^f - 840 {P^f (P^) + 240 q^ (P^) - 1440 q^ {P^f 

- 1260 (Po^)^ - 42672 (P^f {P^f - 11304 q^, {P^f + 5544 (P^i) {P^f - 178416 (Po^)^ 

- 236502 (Po^)^ (Po^)^ - 240408 (Po^)^ (Po^) + 641088 g^ (Po^)^ - 247212 (Po^)^ (Po^)^ 

- 13944 {P^y (P^) - 2520 {P^f (P^) - 15936 {P^Y ql + 996 (Po^)^ + 40068 (Po^)^ + 12 {P^f 

- 165 {P^Y - 2835 (Po^)^ + 84 (Po^)^ - 202986 (Po^)*^ - 2199 {P^Y + ^) 
+^5(577584 go' [P^Y (^o)' " 508032 g' (Poi)^ + 2830464 g' (Poi)^ [P^Y 

- 1 1009376 {PlY {Pof Qo " 3448368 go^ (Po^)^ (P^) + 11760 g' {P^Y (^o ) 
+ 1878912 g^ {P^Y (Pof + 138915 (Po^)"^ (Po^)' + 3873408 {P^Y (PoY Qo 
+ 1 12896 go' {P^Y (^o ) " 780864 g« {P^Y (PoY + 8285760 go^ {P^Y (PoY 

- 47040 go^ (Po^) (P^) - 3321696 g« (P^Y (^o ) " 4235616 (Po') g« (P.'Y 

+ 11690784 (Poi)6 g^ (P')' - 2555280 (Pq^)^ {P^Y 4 " 2006928 g^ {P^Y (^o )' 
+ 1561728 (P2)« go' (Po^) + 9767520 (Po') go^ {P^Y " 19885824 {P^Y (Pof 

- 4906944 gS (Po^)^ {P^Y " 1916880 go^ {P^Y (^o )' " 338688 go' (Po^)^ (Po') 
+ 411264 (Po^)' (Po')' go^ - 6339648 g« (P^Y (PoY + 181440 go' (Po')^ (Po^) 

- 622944 (Po^)^ (Po')^ go' + 5608512 go^ (Pq')' {P,'Y - 26925696 (P^) {P,'Y qI 

- 1750896 go' (Po^)^ (Po')^ + 703080 go^ (Pj)^ (P^) + 16039296 (Po') go^ (Po^)^ 
+ 1300320 g^ (Po^)' (Po')' - 773136 (Po')^ {P^Y - 2721600 (Po') {P^Y 

- 11760 go' (Po')' (Po^) - 570528 g^° (Po^) (Po')^ - 8170848 (Po') go^° (Po^)^ 
+ 4898880 go' (Po')^ (Po') + 2779392 g^ (Po')^ (Po') - 1202544 g^ (Po^)^ 

+ 5450592 go^ (Po')^ (Po')' + 107184 ql {P,') {P^Y + 178752 go' (Po')^ (Po')' 
+ 8013600 g^ (Po')' (Po')' + 9508128 g^ (Po')' (Po')' - 103488 g^ (Po')' (Po')' 

- 7560 go^ (Po') (Po')' - 157248 g« (Po')' (Po') - 1912512 go'° (Po')' (Po')' 
+ 2036160 ql (P^) (Po')' + 5197920 go^ (Pj)' (Po')^ + 41472 go^ (Po')^ 

+ 4110624 go^ (Po')^ (Po')' - 665280 go' (Po')^ (Po') + 278208 go^ (Po')' (Po') 
+ 364896 go^ (Po')^ (Po') - 1014720 go^ (Po')' (Po') - 2610720 go^ (Po')' (Po')' 

- 2550240 ql (Po')' - 9748032 go'° (Po')' (Po')' + 170640 g^ (Po')^ 

- 1201200 go^ (Po')' (Pj)' + 7336224 (Pj)^ (Po')' go' - 1330560 go' (Po')' (Po')' 

+ 53361 (Po')' (Po')' + 6365856 (Po')' (Po')' g^ + 178920 g^ (Po')' + 23520 g^ (Po')' 
+ 120687 (Po')' (Po')' - 840 go' (Po') - 264480 g^ (Po')' + 10976 g^ (Po')' 



37 



- 1496880 {P^Y (P2)4 + 142884 {P^f (P^) + 105 (P^) (pj) + 153120 {P^f 

- 4830 {P^y {P^y + 281472 {P^f ql + 246288 {P^f [P]) + 29568 ql {P^f - 1680 (Pq^)^ [P]) 
+ 360 (Po ) - 3360 ql {P^f - 3360 q^ {P^f - 341376 {P^f {P^f - 56520 (^o )' 

+ 19404 (Po^) (P2)5 - 1784160 (Po^)^ g^ - 1892016 (Po^)^ (Po^)^ - 1923264 (Po^)^ (Po^) 

476448 

+ 5769792 ql {P^f - 1977696 (Po^)^ {P^f - 25982208 (Po^)^ go' - go'° (Po^)^ 



+ 14270256 ql {P^f + 3342528 g^ (Po^)^ + 10644480 g^ {P^f + 9743328 (Po^) {P^f 



111552 (P2)^ (Pq^) - 5040 {P^f (P^) - 143424 (P^Y ql + 4224276 (Po^)^ (P^ 
1239072 



- 18432 g^ {P^Y + go^° (Po')^ + 10877328 {P^Y (^0 )' + 627792 g^ {P^Y 



371 65968 

+ 10724238 {P^Y (Pof + " (^0)' (^0)' + 9821196 (Po^)^ {P^Y + 2083536 (Po^)=^ {P^Y 



- 604608 g^ (Po^)^ + 853230 (Po^)^ (Po^)' + 3486 {P^Y + 140238 (Po^)^ + 15 (Po^)' - 330 {P^Y 

- 5670 (Po^)'^ + 105 (Po^)' - 1623888 (P^Y " 17592 {P^Y + (PoY' + - 

Fso{A){M'o - gMl mo - gmo) = 

3 1n(2) + iln(go^^)-iln((PoY(PoY) 
+g{3 - 32 go' (Po^) + 6 (Po^)' + 15 {P^Y + 24 (P^) (P,')) 

+/(12(Po^)' 

345 

- 240 (Po')^ (Po^) - 45 {P^Y - 32 go^ (Po^)' - ^ (Po')^ + 384 go' (Po^)' (P^) + 96 go^ (Po')' 

- 192 go^ (Po^) (Po') + 192 go' (Po^)^ + ^ " 276 [P^Y {P'of " 144 (Po^)^ (Po') - 96 go' (Po^) 
+ 48 (Po') (Po^) + 30 (Po')') 

)1\5 , 1 non „2 / d1n2 / d2n , /i on ^ d1n4 / d2n2 



y (1 - 2880 go^ {P^Y + 1920 go^ (Po^^ (Po^) + 4212 {P^Y (PoY " 1152 go^ (Pj) (Po^) 

, 



6912 go' (Po^)^ (Po') + 3456 (Po^)' (Po')' go^ + 384 go^ (Po')^ (Po^) + 2304 go^ (Po^)' (Po'^ 



- 3840 go' (Po')^ (Po^)' - 8832 go' (Po^)^ (Po')' - 4224 go^ (Po^) (Po')' - 1728 go^ (Po')^ 
+ 1920 go^ {P^Y (Po) + 6540 (Po')^ (Po^)' + 1920 go^ {P^Y + 960 go' (Po^)^ + 6952 {P^Y (Po? 

- 192 go' (Po^) - 192 go^ (Po^)' + 2160 (Po^)^ (Po') + 72 (Po') (Po^) - 1104 (Po^)' (Po')' 

- 960 (Po')3 (Po^) + 576 go^ (Po')' - 256 q', {P^Y + 5520 (Po^) (Po')^ - 576 (P^Y (^0') + 3485 (P, 
+ 606 (P^Y + 45 (Po')' - 690 (Po')^ - 180 (Po^)^ + 18 (Po^)' + 1728 go^ (Po^)^) 

+^7^(333696 go' (Pj)^ (Po')=^ - 20160 go' (Po^)^ + 269568 go' (Po^)^ (Po')' - 82944 go^ (P^Y (^0) 
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+ 5760 ql {P^f (P^) + 25272 {P^Y {P^f - 4032 q', (P^) (P^) - 110592 {P^Y {P^f q', 

- 27648 q', {P^f (P^f - 43008 q', (P^)' (Pof " 48384 q', (P^)' (P^) + 27648 {P^f {P^f q', 
+ 3072 q', {P^f (P^) + 209280 q^ {P^f {P^Y + 20736 q', {P^f (P^) + 80640 g« {P^f {P^f 

- 246600 (P^f (P^f - 100224 (P^) (p^i)5 ^4 ^ ^73800 q^ (P^f (P^) + 2304 g« (Pq^)^ 

+ 84864 g6 (p^i) (p2)4 _ 2688O (^o')' (^0 " 61824 ql (P,')' (P^)' - 38016 (^o') (Pof 

+ 34560 g« {P^f (P^) - 28416 g« (P^) {P^f " 13824 go^ (Po^)^ + 15360 go^ (Po^)^ (Po^) 

+ 11136 go^ (Po^)^ (Po^) - 99840 go^ (P^)' iPo')' " 35712 ql {P^f + 6912 g« {P.^ 

+ 67584 g^ (Po')^ (Po^)" + 88320 g^ {P^f {P^f + 39240 (Pq')^ (Po')' + 17280 g^ (Po^)^ 

+ 2880 go^ {P^f + 41712 (P^f (P^f - 320 go^ (Po^) + 45312 go^ (P^)' - 672 go^ {P^f 

- 238374 (Po^)^ (Po^)^ + 12960 (P^)' (P^) + 96 (P^) (Po^) + 9984 gg^ {P^f - 2760 (P^^ {P^f 

- 2400 (Po^)^ (Po^) + 2016 ql {P^f - 214260 {P^f {P^f - 2304 g^ {P^f + 33120 {P^) {P^f 

- 61056 (Po^)^ g^ - 163224 (Po^)^ (Po')=^ - 43632 {P^Y {P^) + 58176 go' (Po')' + 60 (Po')' 

- 94500 (Po^)*' {P^f - 167280 {P^^ {PD - 1440 (Po^)^ (Po^) + 20910 (Po^)^ + 3636 {P^f 

- 1725 (Po')^ - 450 (Po^)^ + 24 (Po^)' - 10665 {P]f - 91140 (Po')*^ + ^ + 13824 g^ {P^Y) 
+/(3003264g2 (P^f (P^f - 80640 g' {P^f + 2426112 g' {P^f {P^f 

- 13057920 (Po^)'^ (P2)3 ql - 912384 go^ {P^^ (Po) + 13440 g' (Po^)' (P^) 

- 2666496 q^ {P^f {P^f + 88452 (Po^)^ {P^f - 2676480 (Po^)' (Po')' gg' 

- 2239488 g« {P^f {P^f + 4918272 go^ (Po^)^ {P^^ - 10752 go^ (Po^) (Po') 
+ 476160 g« (P^)' (P^) - 1244160 (P^) ql (P^f + 5985792 {P^f q^ {P^f 

- 1105920 (Po^)^ (Po')' go^ - 304128 go^ {P^f {P^f + 2059776 (Po') go^ (Po^)^ 

- 9072000 (Po^)' (Po')' ql - 2276352 g^ (Po^)^ (Po')' - 430080 ql {P^f {P^f 

- 193536 go' (Po^)^ (Po') + 124416 (Po^)' (Po')' go^ - 602112 g« (Po^)' (Po')^ 
+ 13824 go^ (Po')^ (Po^) - 15255936 (Po^)^ (Po')^ go' + 2366208 go^ (Po')' (Po^)^ 

- 4886784 (Po') (Po^)« go' + 1883520 go' (Po^)^ (Po')' + 103680 go^ (Po^)' (Po') 
+ 4745088 (Po') go (Pq )' + 967680 g^ (Po^)' (Po')' - 1972800 (Po')^ {P^f 

- 1002240 (Po') (Po^)^ go' - 723456 g^° (Po^) (Po')' - 958464 (Po') g^° (Po^)' 

+ 1555200 go' (Po^)^ (Po') - 836736 go' (Po')' (Po^) + 27648 g« (Po^)' + 663552 go' (Po')^ (Po^)' 
+ 933504 go^ (Po^) (Po')' - 107520 go' (Po')=^ (Po^)' + 2390016 go' (Po^)^ (Po')^ 
+ 6532608 go' (Po^)^ (Po')^ - 247296 go' (Po^)^ (Po')' - 190080 go^ (Po^) (Po')' 
+ 414720 g^ (Po^)^ (Po') + 233472 g^° {P^f {P^f - 340992 g^ (Po^) (Po')=^ 
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- 1826304 ql {P^f {P^Y - 62208 q', {P^Y + 567552 q', {P^Y {P^f + 69120 q', (P^)' (^o ) 

- 2517120 q', (P^f iP,') + 111360 q', (P^)' (P,') - 998400 q^ {P^f {P^f 

- 392832 q^ (Pj)^ + 1812480 {P^f {P^f + 82944 (P^)^ + 743424 q^ {P^f {P^f 

- 11836800 {P^y (P2)5 g2 + 794880 q^ {P^f {P^f + 137340 (P^)^ {P^f 

- 6856320 (Po^)3 (P^)^ q^ + 86400 (Poi)=^ + 6720 ql {P^f + 145992 {P^f {P^f 

- 480 ql (Po') + 453120 {P^f - 1792 g^ (Pq^)' - 1906992 (P^)^ (P^)^ + 45360 {P^f (P^) 

+ 120 (P2) (Pgi) + 109824 g^ (P^)^ - 5520 {P^f {P^f + 5832960 {P^f {P^) - 4800 {P^f {P^) 
+ 5376 g^ {P^f - 1714080 (Po^)^ (Po^)^ - 11520 g^ (P^)^ + 115920 {P^) (P^)^ 

- 610560 {P^f q^ - 1305792 (Pq^)^ (P^)^ - 349056 (Pq^)^ (Pq^) + 523584 g^ {P^Y 

- 756000 {P^r {P^r - 1365120 (Po^)^ ql - g^° (Pq^)^ + 2039040 q', {P^f 

- 622080 g^ {P^ f + 374400 g^ {P^Y + 1023840 (Pq^) {P^Y - 1338240 (Pq^)^ (Pq^) 

- 2880 {P^Y (Pq) + 10791720 {P^Y (PoT - 395520 g^ {P^Y - 67584 g^^^ {P^Y 

49990752 

+ 4617216 {P^Y (Po? - 1434624 go^ {P^Y + 2490696 {P^Y {Pof + % (^o)' (n')' 

+ 7313976 {P^Y (Pq)^ + 9979200 (Po^)^ {P^Y - 311808 g^ {P^Y + 8160000 {P^Y (Pq ? 

+ 73185 {P^Y + 12726 (Pq^)^ + 75 {P^Y - 3450 (Pq^)^ - 900 {P^Y + 30 {PqY - 85320 (Pq^)*^ 

- 729120 {P^Y + 2769360 (Po')^" + ^ + 12^^ (p^i)io + 62208 q^ {P^Y)- (H.2) 

I Numerical values of extrema 

The positions of the extrema having the lowest free energy are shown in Table |^. Here Pq, 
Pq and go are related to Mq, Mq and mo as in appendix ^ Since we consider lower orders 
here, we can not expect clear plateaus. However at lower orders we expect that the extremum 
with the lowest free energy lies in the plateau and approximates the free energy well. We can 
check it is indeed the case for the 0^ toy model. 
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ansatz 
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1st ordpr 


3rd ordpr 


5th ordpr 




pi 


n /1 1 nc;c;n7Qnc; 
u.^tiuoou ( yuo 




u.^t^tuuzzooyo 
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go 
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go 
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Table 2: The positions of the extrema for each order. 



[3] J. Maldacena, "The large limit of superconformal field theories and supergravity," 
Adv. Theor. Math. Phys. 2, 231 (1998) [Int. J. Theor. Phys. 38, 1113 (1998)] ||arXiv:hep-| 
th/9711200|. 



[4] O. Aharony, S. S. Gubser, J. Maldacena, H. Ooguri and Y. Oz, "Large N field theories, 
string theory and gravity," Phys. Kept. 323, 183 (2000) ||arXiv:hep-th/9905111 |. 



[5] T. Banks, W. Fischler, S. H. Shenker and L. Susskind, "M theory as a matrix model: A 
conjecture," Phys. Rev. D 55, 5112 (1997) | |arXiv:hep-th/96l0043|| . 

[6] N. Ishibashi, H. Kawai, Y. Kitazawa and A. Tsuchiya, "A large-N reduced model as 
superstring," Nucl. Phys. B 498, 467 (1997) ||arXiv:hep-th/96l2ll5|| . 

[7] L. Smolin, "M theory as a matrix extension of Chern-Simons theory," Nucl. Phys. B 591, 
227 (2000) |^rXiv:hep-th/0002009 |. 

[8] T. Azuma, S. Iso, H. Kawai and Y. Ohwashi, "Supermatrix models," Nucl. Phys. B 610, 
251 (2001) |^rXiv:hep-th/0l02l68 |. 



[9] L. Smolin, "The exceptional Jordan algebra and the matrix string," ||arXiv:hep^ 
th/0l04o50| . 



[10] H. Aoki, S. Iso, H. Kawai, Y. Kitazawa, A. Tsuchiya and T. Tada, "IIB matrix model," 
Prog. Theor. Phys. Suppl. 134 (1999) 47 [|arXiv:hep-th/9908038 |. 



[11] J. Nishimura and F. Sugino, "Dynamical generation of four-dimensional space-time in 
the IIB matrix model," ||arXiv:hep-th/0111102|| . 

[12] D. Kabat and G. Lifschytz, "Approximations for strongly-coupled supersymmetric quan- 
tum mechanics," Nucl. Phys. B 571, 419 (2000) ||arXiv:hep-th/9910001|| ; 



41 



S. Oda and F. Sugino, "Gaussian and mean field approximations for reduced Yang-Mills 
integrals," JHEP 0103 026, (2001) ||arXiv:hep-th/0011175|| ; 

F. Sugino, "Gaussian and mean field approximations for reduced 4D supersymmetric 
Yang-Mills integral," JHEP 0107 014, (2001) ||arXiv:hep-th/0105284|| . 

[13] For example, see R. Fukuda, M. Komachiya, S. Yokojima, Y. Suzuki, K. Okumura and 
T. Inagaki, "Novel use of Legendre transformation in field theory and many particle 
systems," Prog. Theor. Phys. Suppl. 121 (1995) 1. 

[14] E. Brezin, C. Itzykson, G. Parisi and J. B. Zuber, "Planar Diagrams," Commun. Math. 
Phys. 59, 35 (1978). 

[15] J. Nishimura and G. Vernizzi, "Spontaneous breakdown of Lorentz invariance in JIB ma- 
trix model," JHEP 0004 015, (2000) [arXiv: |hep-th/000322l ; 

"Brane world from JIB matrices," Phys. Rev. Lett. 85, 4664 (2000) [arXiv: [hep- 
th/0007022| ; 

J. Ambjorn, K.N. Anagnostopoulos, W. Bietenholz, T. Hotta and J. Nishimura, "Monte 
Carlo studies of the JIB matrix model at large N," JHEP 0007 Oil, (2000) [arXiv: |Eep= 
th/0005T47i ; 

K.N. Anagnostopoulos and J. Nishimura, "A new solution of the complex action problem 
and the dynamical space-time in the HB matrix model," [arXiv: [hep-th/0108041|| ; 
Jun Nishimura, "Exactly Solvable Matrix Models for the Dynamical Generation of Space- 
Time in Superstring Theory," Phys. Rev. D 65, 105012 (2002) [aeXiv: |hep-th/ 01 080701 . 



42 



